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Abstract 

The moderately large measured value of #13 signals a departure from the approximate two-flavour 
oscillation framework. As a consequence, the relation between the value of #23 i n nature, and the 
mixing angle measured in Vn disappearance experiments is non-trivial. In this paper, we calculate 
this relation analytically. We also derive the correct conversion between degenerate values of #23 
in the two octants. Through simulations of a v„ disappearance experiment, we show that there are 
observable consequences of not using the correct relation in calculating oscillation probabilities. 
These include a wrong best-fit value for 023) and spurious sensitivity to the octant of #23. 
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I. INTRODUCTION 



Neutrino oscillation physics has entered an era of precision measurements. Since 2011, 
the reactor experiments Daya Bay [1], Double Chooz [2, 3] and RENO [4] and superbeam 
experiments MINOS [5] and T2K [6] have measured a non-zero value of 9 13 [7-9]. Analyses 
of world neutrino data [10-12] have given us the value sin 2 2#i 3 ~ 0.1, which is moderately 
large. Daya Bay has measured sin 2 2# 13 = 0.089 ± 0.011 [13], which is the most precise 
measurement till date. These measurements have established that #13 is non-zero at more 
than 5cr confidence level. 

The solar mixing angle # 12 and mass-squared difference A 2 i (A^- = m 2 — m 2 ) have been 
measured very accurately by SNO [14] and KamLAND [15], respectively. Their current 
best-fit values are sin 2 12 = 0.32 ± 0.017 and A 2 i = (7.62 ± 0.2) x 10~ 5 eV 2 [10]. MINOS 
[5] has measured the mixing angle 6^ and mass-squared difference A pfl with a precision 
of a few percent. We use the subscript fifi to indicate that these parameters are measured 
from observations of muon neutrino disappearance. The values of these parameters from 
MINOS are sm 2 26 IHi > 0.90 (90% C.L.) and \A^\ = (2.32±0.1) x 10~ 3 eV 2 [16]. In the two- 
flavour oscillation scenario (neglecting the small parameters A 2 i and #13), the muon neutrino 
survival probability P MM depends on sin 2 29^ = sin 2 2# 23 and |A MM | = |A 31 |. Because this 
probability depends on the magnitude but not the sign of |A 31 |, we cannot determine the 
neutrino mass ordering or hierarchy. Moreover, this function is symmetric about # 2 3 = 45°, 
which gives rise to the octant degeneracy. Thus, the currently unknown parameters in 
standard three-flavour neutrino oscillation physics are - (a) the neutrino mass hierarchy 
(normal hierarchy (NH): A w > or inverted hierarchy (IH): A w < 0), (b) the CP- violating 
phase 5cp and (c) the octant of 9^ (lower octant (LO): 9^ < 45° or higher octant (HO): 
&w > 45°). 

It is possible for the mass hierarchy to be determined by the upcoming experiment NOz/A 
itself, if the value of 5c p in nature is in the favourable range [17]. Combined data from 
multiple experiments [18-20], experiments with longer baselines [21, 22] and atmospheric 
neutrino experiments [23, 24] can also determine the hierarchy. The measurement of 5c p is 
difficult but possible, thanks to the non-zero value of #13. This requires very intense beams 
at short baselines [25]. In this work, we concentrate on the precision measurement of the 
atmospheric mixing angle. 
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#23 is the largest of the mixing angles in the leptonic and quark sectors. Its near-maximal 
value is indicative of a symmetry of nature in the fi — r sector [26, 27]. In many theoretical 
models, the deviation of 6*23 from maximality is related to the deviation of 6*13 from zero [27]. 
Thus, the precision measurement of this angle, and determination of its octant can play an 
important role in constructing new models of physics. 

In this paper, we discuss the measurement of #23 and its octant, particularly in light of 
moderately large #13, which gives rise to three- flavour effects. In Ref. [28, 29], the authors 
had shown that due to three-flavour effects, the mass-squared difference A M/J measured in 
muon disappearance experiments is not A 31 , but a linear combination of A 31 and A 2 i- The 
same calculation also indicates that the mixing angle 6^ measured in these experiments is 
not the same as #23- While this result has been seen in the literature before (in Ref. [30]) and 
more recently in Ref. [12], we present a detailed study of this effect in this paper. We pay 
particular attention to the effect of choosing the 'wrong' definition (# 23 = # wt ) in analyses. 
In Section II, we have outlined the calculation that indicates the relation between A MAt , 9^ 
and A3i,#23- We have also found analytic expressions relating deviations from maximality 
in the two octants. In this work, our aim is to highlight a physics point, rather than study 
the capability of any particular experiment. However, in order to make our point clearer, 
we have presented the results of some simulations, in Section III. Finally, in Section IV, we 
have summarized our findings. 

II. CALCULATIONS 

The P w oscillation probability in the three-flavour scenario (ignoring matter effects) is 
given by 

p m =\ - 4i^ 3 ri^ t ii 2sin2A 3i 

- 4|U" Ai3 | 2 |^2| 2 sin 2 A 32 

- 4|^ 2 | 2 |^ 1 | 2 sin 2 A 21 , (1) 

where we have used the shorthand notation Ay = AijL/AE. Here, U a i are the elements of 
the 3x3 PMNS matrix. This probability is sensitive to all six standard oscillation parameters 
- three mixing angles, two mass-squared differences and the CP phase. In interpreting the 
result of oscillation data in terms of two-flavour oscillations, we attempt to express the 
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probability in the simple form 



P w = 1 - sin 20 w sin A w 



(2) 



that involves only two parameters. Recasting the full six-parameter formula as a simple 
two-parameter formula results in a non-trivial relation between A w and A31, and between 
and #23- I n this calculation, we will consistently retain only terms upto linear order in 
A21 (since A21 << A31). Ignoring the last term in Eq. 1, we have 



l-P, 



/</< 
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|t/ M 3| 2 |^i| 2 sin 2 A 31 + |^3| 2 |^2| 2 sin 2 A 
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Here we introduce the notation 
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Note that U\ + U2 — 1. This notation lets us write 
l-P, 





U, 2 


to 




2 + 


Uu2 





MM 



U^\ 2 (\U^\ 2 + If/^l 2 ) Ui sin 2 A 3 i + u 2 sin 2 A 32 



Using the fact that A 32 = A31 — A 2 i, a little algebra gives us 

1 - cos 2 A 3 i - 2w 2 A 2 i sin 2 A 



l -^ = \\u^\ 2 {\u^ + \u^) 



Ml 



We rewrite this as 



l-P, 



mm 



^rocvr+iM 2 ) 



cos /3 cos 2 A31 + sin j3 sin 2A 3i 



1 + 4u 2 A 2 1 



where 



cos f3 



^ 3 | 2 (|^i| 2 + |^ 2 | 2 ) 



1 - cos(2A 31 -P)a/1 + 4m 2 A 



- 2 ^ii 



1 + 4n 2 A 2 x 



sin (3 



2w 9 A 



2^21 



1 + 4n 2 A 2 x 



(3) 



(4) 



(5) 



(6) 



We ignore the quadratic term in the square root, and we note that j3 = tan 1 (2m 2 A 2 i) ~ 
2m 2 A 2 i. This gives us our final result 

P m = 1 - A\U^\ 2 (1 - |^ 3 | 2 ) sin 2 (A 3 i - n 2 A 2 i) . (7) 
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On comparing with the two-flavour formula, we can make the association 



A 3 i = A w + u 2 A 2 i , 
which is the result expressed in Ref. [28]. Moreover, we also find that 



sin 2 20 w = 4|C/ M3 | 2 (1-|[U 2 ) 

= 4 cos 2 0i3 sin 2 23 (1 - cos 2 6 13 sin 2 23 ) ■ 

This means that 

sin 023 = sin MM / cos 0i3 or sin0 2 3 = sin(90° - w )/cos0 13 . (8) 

In other words, given a value of sin 2 20 w , there are two degenerate allowed values of 8 m : 
6^° (in the lower octant) and 0^° (in the higher octant). These are related by 



C = 900 - C • ( 9 ) 



Then, the corresponding values of 023 are given by 

Tn sin Off „ n sin0^° , 
sin0 2 L 3° = £L ; sin0f 3 ° = ^- . 10 

16 COS 013 2i COS 013 V ' 

Thus, if 0i3 = 0, we have simply Q^z = an d ^23° = ^au '■ However, for non-zero 0i 3 , the 
relation between Q\° and 9^ 3 ° depends on the value of i3 . 

We have shown this feature in Fig. 1. Along the x-axis, we have different values of sin 2 29^ 
as measured using the P MM channel. For a particular value of 0i 3 , each MM corresponds to two 
values of 023- These two values have been plotted along the y-axis. Thus given a value of 0i3 
we have a curve that maps one value of MM to two values of 023- To show the effect of 0i 3 , we 
have varied it in the range sin 2 20i3 G [0.0, 0.2]. The solid black curve is for sin 2 20i3 = 0.1. 
The spread in i3 gives the shaded band in the figure. 

The main point to be noted here is the asymmetry of the band about the 023 = 45° 
line. When sin 2 20i 3 = 0, we see that the curve is symmetric about 023 = 45°. However, 
for larger values of 0i 3 such as sin 2 20i3 = 0.2, we lose this symmetry. For instance, if 
sin 2 20i3 = 0.2 and if the MINOS measured value of sin 2 20 MM = 0.98, then the values of 023 
in the two octants are 40.5° and 52.7°. These two angles are not complementary. In other 
words, it is MAt that goes to 90 — MM under the octant degeneracy; but 02 3 ° and 0^° are 
not complementary angles - the exact change depends on the value of 0i3. Henceforth, for 



w 
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CD 

D) 
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0.92 0.93 0.94 0.95 0.96 0.97 0.98 0.99 



sin 2 29. 



FIG. 1: Allowed values for #23 (y-axis) in degrees, for a given measurement of sin 2 29^^ (x-axis). The 
dashed (blue) and dotted (green) curves are for the extreme cases sin 2 2#i 3 = 0.2 and 0.0, respectively. The 
solid (black) curve is for sin 2 29 1$ =0.1. 



convenience, we will refer to 0^ and 0^ corresponding to a given measurement of sin 20 w 
as being /i-complementary . 

As an interesting aside, we note that (given sin 2 20 13 = 0.1) if 44.3° < MM < 45.7°, i.e. if 
sin 2 29 > 0.9993, then both allowed values of #23 are greater than 45°. In such a case, there 
is no ambiguity in the octant of #23 - it is necessarily in the higher octant. However, current 
experiments do not have the precision to distinguish sin 2 26^ = 0.9993 from sin 2 20 w = 1. 



Therefore, this point is purely of academic interest. 

Hereafter, we assume that sin 2 20 13 = 0.1. When sin 2 29 



we get # 23 



, -= 1.0, i.e. 9% = 9™ = 45°, 
- 0^° = 45.75° from Eq. 10. This is also seen from Fig. 1. For a given value 



of sin 2 20 MA1 7^ 1.0, the two allowed values of 9^ have equal and opposite deviations from 
45°, i.e. 9^9 < 45° and 9^9 > 45°. However, the two /i-complementary values of 6*23 lie on 
opposite sides of 45.75°. Since 

sin 6j£ = sin 0^° cos 13 and 

sin = cos 9 L £ = sin 0g° cos0 13 , 
one can eliminate 9z9 between these two equations. This gives us 



sin 9g° 



1 



COS 2 013 



sin 2 9h? 



23 > 



(11) 



which is a handy equation to switch from one octant to another. For example, if #23 = 40°, 
this equation tells us that the corresponding value of 9^° is 51.54°. 

We can try to recast Eq. 1 1 in terms of deviations from maximality, rather than in terms 
of the angles themselves. To this end, we define 

S^° /HO = 9 2 , LO ' HO - 45° . 

Assuming small deviations, we can linearize Eq. 11. This gives us 

8g° = sin 2 9 13 - = 1.43° - 4? , (12) 

which implies 

9g° = 91.43 - e£ . (13) 

We will use these relations to interpret the results of our simulations. 

In Ref. [30], the authors have argued that a precise measurement of 6*23 is difficult in the 
vicinity of 45° because sin 2 2^23 attains a maximum here, making A(sin 2 2#23) very small. 
Therefore it is worth discussing whether this ^-effect can have experimentally observable 
consequences. In order to directly observe a shift of 0.75° at 9^ = 45° (#23 = 45.75°), we need 
a precision of 0.001 in our measurement of sin 2 26*23. This is far beyond our experimental 
reach. However, if = 40°{9^ = 40.63°), the precision in sin 2#23 required to observe this 
shift is around 0.007, which may be achievable at future facilites. Moreover, this effect can 
be felt indirectly, as an artificially enhanced/reduced sensitivity in oscillation experiments. 
We illustrate this through simulations, in the next section. 

III. SIMULATIONS 

As we have mentioned before, our aim is to illustrate the difference between choosing the 
'wrong' definition: #23 = 9^ and the 'right' definition: #23 = sin _1 (sin 9^/ cos #13). In order 
to show the experimentally observable effects of this choice, we have simulated the NOz/A 
experiment using the GLoBES package [17, 31-36]. In this section, we discuss the results of 
our simulations. 

We have used the standard NOz/A setup [17], with a 14 kton totally active scintillator 
detector placed 810 km away from the NuMI beam at a 14 mrad off-axis location. The 
beam, with a power of 0.7 MW, is assumed to run for three years each with neutrinos and 
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FIG. 2: i^, disappearance event rates for 6*23° = 40° are shown as a solid (red) histogram in both panels. 
Superimposed on this, as a hatched (green) histogram are the event rates for the corresponding #|| using 
the 'wrong' definition (left panel) and 'right' definition (right panel). For the 'right' definition, the event 
rates show a good match. 

antineutrinos. We have taken the energy resolution for v yi to be 0. 06 ^/E(GeV) [36]. Back- 
grounds from NC events have also been taken into account. For the oscillation parameters, 
we have chosen sin 2 12 = 0.304, sin 2 26 l3 = 0.1, A 21 = 7.6 x 10^ 5 eV 2 , A w = 2.4 x 10~ 3 eV 2 
and 5cp = 0, unless specified otherwise. 

In Fig. 2, we have plotted the event rates from the muon disappearance channel. In the 



LO 



left panel, we have chosen the 'wrong' definition of 623, so that 9 23 



40° implies 9, 



HO 
23 



50°. 



This gives us a difference in the number of events. However, on using the 'right' definition 
(^23° = 51.54°), we find that the event rates match, as seen in the right panel. 

Having showed the difference due to our choice of #23 at the event level, we now proceed to 
do so at the level of x 2 . Throughout our simulations, we have used the values specified above 
as the true values of oscillation parameters. We have varied the test values of the parameters 
in the following ranges: sin 2 26 13 e [0.07, 0.13], W G [35°, 55°], \A^\ e [2.05, 2.75] x 10" 3 eV 2 
and 5qp £ [0,27r). The solar parameters A21 and 612 have been kept fixed in this analysis, 
since the effect of their variation is small. The mass hierarchy is assumed to be normal. 
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23 (test) [higher octant] 
55 54 53 52 51 50 49 48 47 46 




36 37 38 39 40 41 42 43 44 45 
23 (test) [lower octant] 



23 (test) [higher octant] 
55 54 53 52 51 50 49 48 47 46 




36 37 38 39 40 41 42 43 44 45 
23 (test) [lower octant] 



FIG. 3: x 2 f° r sensitivity to the value of #23 from muon disappearance from NO^A . True = 40°. The 
solid (pink) curve gives \ 2 for test 9^ in the lower octant (read the lower x-axis), wile the dotted (green) 
curve gives % 2 for test 6^ in the higher octant (read the upper x-axis). In the left panel, we have chosen 
the 'wrong' definition of $23 while in the right panel we have chosen the 'right' definition of 6*23. 

In Fig. 3, we have plotted the sensitivity to #23 for true 9^ = 40°. As the test value of 
9^9 increases in the lower (true) octant from 35° to 45°, 9%f increases from 35.53° to 45.75°. 
This range is shown on the lower x-axis. Correspondingly, 9^° decreases from 55.60° to 
45.75°. This range, for the higher (false or degenerate) octant, is shown on the upper x- 
axis. The advantage of using these double-axes is that values of #23 along a vertical line are 
/i-complementary . We have used the linearized relation in Eq. 13 to plot these axes. For 
values of test 9^ in the lower octant, y 2 has been plotted as the solid (pink) curve. The 
lower x-axis should be used to read the #23 values for this curve. For values of test 9^ in 
the higher octant, x 2 has been plotted as the dotted (green) curve. The upper x-axis should 
be used to read the #23 values for this curve. 

The following features are visible in Fig. 3: (a) For 9^ = 40°, the best-fit point is seen at 
the corresponding value of #23 which is 40.63° (b) If we use the 'wrong' definition of #23, the 
minima in the two octants do not appear for /i-complementary values of #23- This is seen in 
the left panel. In the right panel, we have used the 'right' definition. As a result, we find 
that the minima coincide, (c) The sensitivities in the two octants coincide at #23 = 45.75°, 
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23 (test) [higher octant] 
55 54 53 52 51 50 49 48 47 46 



23 (test) [higher octant] 
55 54 53 52 51 50 49 48 47 46 





36 37 38 39 40 41 42 43 44 45 
23 (test) [lower octant] 



36 37 38 39 40 41 42 43 44 45 
23 (test) [lower octant] 



FIG. 4: x 2 f° r sensitivity to the value of #23 from muon disappearance from NO^A . True 



45°. The 



solid (pink) curve gives \ 2 for test 9^ in the lower octant (read the lower x-axis), wile the dotted (green) 
curve gives % 2 for test 6^ in the higher octant (read the upper x-axis). In the left panel, we have chosen 
the 'wrong' definition of 023 while in the right panel we have chosen the 'right' definition of 6*23. 

rather than 823 = 45°. This is expected from the analytic calculations presented in the 
previous section. 

Figure 4 is similar to Fig. 3, but with true 8^ = 45°. Note that compared to the 'right' 
definition x 2 , the 'wrong' definition \ 2 is more in the lower octant, and less in the higher 
octant. Clearly, using the 'wrong' definition gives us an incorrect value of \ 2 ■ Thus, if x 2 
from the 'wrong' definition is used to set a prior on 6*23 for future experiments, it will give a 
spurious indication of the confidence with which certain 6*23 values are allowed. This also has 
implications for octant sensitivity studies. Since 823 and 90° — 823 are not /i-complementary , 
X 2 (@23', 90° — 823) does not give the correct octant sensitivity. Therefore, care must be taken 
to use the 'right' definition in simulations. 

Finally, in Fig. 5, we have plotted the allowed values of 823 for all possible true values of 
#23 in the range [36°, 55°]. For each true value, the allowed test values can be read from its 
corresponding vertical line. The 68%, 90% and 95% C.L. contours are also shown. For a 
given true value of #23, we find an allowed range of test 823 in each octant. For the 'right' 
definition, we see that the false degenerate solution lies along the line 8^° = 91.43° — 8^ 
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36 38 40 42 44 46 48 50 52 54 36 38 40 42 44 46 48 50 52 54 

e 23 (true) e 23 (true) 

FIG. 5: Allowed test values for 6*23 (y-axis) in degrees, for a given true value of 6*23 (x-axis). 68%, 90% 
and 95% C.L. contours are also shown. As before, the left(right) panel is plotted using the 'wrong(right) ' 
definition of #23. With the 'right' definition, we see the expected behaviour - the false degenerate solution 
lies along the line O^f = 91.43° — 6*23° (shown by the straight diagonal (green) line). 

(Eq. 13). But if we choose the 'wrong' definition, then we get degenerate solutions that are 
not /i-complementary , along with incorrect values of x 2 . 

In this paper, we have only presented the results for the case where NH is the true 
hierarchy, and true Sep = 0. However, we have verified that these results hold for IH, and 
for a number of values of 5c p- 



IV. CONCLUSIONS AND SUMMARY 



In this study, we have discussed the effect of three-flavour mixing on our measurement 
of 6*23. We have calculated the relation between the value of #23 in nature, and 9^ that 
is measured in muon disappearance experiments. The difference between these two num- 
bers is significant, in light of the measured value of #13. Using the 'right' definition of #23 
(incorporating the ^13-effect), we have found the allowed values of #23 in the two octants, 
corresponding to a single measured value of sin 2 29^. We know that complementary values 
of 9^ in the two octants are related by 9^® = 90° — 9jj9. But for a given measurement of 
sin 2 20 w , we have found that the corresponding values of 6*23 in the two octants are related 
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by 6^° = 91.43° — #23°. We have called this relation /i-complementarity. The exact form of 
this equation comes from the value of #13. 

Through simulations, we have showed that using the 'wrong' definition gives us degen- 
erate fits at non-/x-complementary values of #23 . Consequently, the muon disappearance 
analysis can give an incorrect prior on #23 • in determining octant sensitivity, if the 'wrong' 
definition is used, one can get an incorrect value of x 2 . Therefore, we advocate the use of 
the 'right' definition (as given in Eq. 10) for calculating the oscillation probability in sim- 
ulations. However, priors from previous experiments should be added in terms of sin 2 2# MM 
- the quantity that is measured in the muon disappearance experiments. Our results and 
conclusions have been found to hold true for both hierarchies and for values of 5c p in its 
entire allowed range. 
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